Electrostatic cluster formation in lipid monolayers 
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We study phase separation in mixed monolayers of neutral and highly negatively charged lipids, 
induced by the addition of divalent positively charged counterions. We find good agreement between 
experiments on mixtures of PIP2 and SOPC and simulations of a simplified model in which only the 
essential electrostatic interactions are retained. Thus, our results support an interpretation of PIP2 
clustering as governed primarily by electrostatic interactions, in which divalent ions such as calcium 
mediate an effective attraction between like-charged lipids. Surprisingly, the mediated attractions 
are strong enough to give nearly complete phase separation, so that clusters can even form when 
the overall concentration of PIP2 is low, as is the case in the cell membrane. 



PACS numbers: 87.14.Cc, 87.15.nr, 87.16.dt, 41.20.-q 

One of the truisms of biology is that electrostatics play 
a relatively unimportant role in determining structure on 
scales above the Debye screening length, which is on the 
nanometer scale. Highly charged biomolecules such as 
DNA pQ and actin [2] form an exception to this rule by 
aggregating into large bundles in the presence of multi- 
valent ions. Another biomolccule that carries a very high 
negative charge density is the membrane lipid PIP2. De- 
spite the fact that its concentration in the membrane is 
extremely low (of order 1%), this lipid plays a critical role 
in many processes involving the cell membrane, including 
cell division [3], exchange of chemicals with the environ- 
ment through endocytosis and exocytosis [I], and cell 
motility [5]. Evidence exists that PIP2 forms clusters [5] 
at the sub-micron scale under roughly physiological con- 
ditions. It has been conjectured that this clustering is 
crucial to its effectiveness at such low overall concentra- 
tion laig. 

Various mechanisms for the clustering have been pro- 
posed, including PiP2-protein interactions [9J [10] , exclu- 
sion from cholesterol-enriched ordered domains [51 ITT] , 
and hydrogen bonds |121 113) . However, recent experi- 
ments showed that PiP2-clusters are induced simply by 
adding calcium or other divalent ions [6l [14]. This 
raises the question of whether a purely electrostatic, 
counterion-mediated mechanism is consistent with exper- 
imental observations of calcium-induced PiP2-clustering. 

In this letter, we study ion-induced clustering of neg- 
atively charged lipids as a function of lipid charge and 
ion properties, both numerically and experimentally We 
conduct simulations on a model designed to retain only 
the most critical features of the electrostatics and com- 
pare the results to experiments on Langmuir monolayers 
of a mixture of PIP2 with neutral lipids with added diva- 
lent salts. We find semiquantitative agreement between 



the simulations and experiments, with similar trends for 
the dependence on lipid charge and ion size, suggesting 
that multivalent-ion-mcdiatcd attractions arc indeed re- 
sponsible for the observed clustering. 

Counterion-mediated attractions are the collective re- 
sult of a near-cancellation of repulsive and attractive in- 
teractions between like and unlike charges, respectively, 
in strongly correlated electrostatic systems. They are not 
captured in mean-field theory |15j and are typically quite 
weak [16 . Here, we show that counterion-mediated at- 
tractions are surprisingly strong so that phase separation 
is nearly complete even at reasonable values of the PIP2 
charge, implying that clustering can occur even at very 
low PIP2 concentration. 

Experiments — We look for phase separation using vi- 
sual analysis of fluorescence micrographs of mixed lipid 
monolayers prepared in a Langmuir trough (Kibron. Inc.) 
and imaged on an inverted epifluorescence microscope. 
We use a molar fraction </>pjp of L-a-phosphatidylinositol- 
4,5-bisphosphate (PIP2, Avanti) in a monolayer other- 
wise consisting of l-stearoyl-2-oleoyl phosphatidylcholine 
(SOPC, Avanti). SOPC is zwitterionic, and known to be 
net-neutral over a wide range of pH. Part of the PIP2 
(0.5 mol% of the total lipid content) is replaced by a flu- 
orescently labeled analogue (Bodipy FL-pip 2 , Echelon, 
Inc.). A lipid monolayer is formed on a buffered sub- 
phase (10 mM hepes, 100 fjM edta, 5 mM DTt) by addi- 
tion of the lipids dissolved in a 2:1 chloroform/methanol 
mixture to the air-water interface. After formation, the 
monolayer is visualized by epifluorescence to determine a 
baseline frequency of structural inhomogeneity. Divalent 
salts CaCl2 or MgCl2 were added at 1 mM to the sub- 
phase, followed by gentle mixing to avoid disrupting the 
monolayer. We image the monolayers about 50 minutes 
after cation addition to allow sufficient time for domain 
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FIG. 1: Snapshots of the experiment (a) and simulation (b) at <f>pip = 25%. (a) The (inverted) epifluorescence micrograph, 
taken 50 minutes after mixing ImM CaCl2 into the subphase (pH 7.4). PlP2-rich domains are shown as very dark green spots, 
(b) The simulation (piP2-charge qpip = —4, divalent ion radius Rci = 2 A) after 3.5 ps of coarsening. Charged and neutral 
lipids are drawn green and light grey, respectively, and divalent ions that are close to the lipid monolayer are plotted in dark 
red. (c) Strength (shaded contours) and direction (streamlines) of the electric field around a string-like domain (geometry taken 
from the simulation), illustrating that further growth of the domain is likely to occur at the end. 



coarsening. The existence of bright spots at a higher fre- 
quency than baseline serves to determine the existence of 
PiP2-rich clusters. 

We perform this procedure for a range of </>pip-values 
and several pH values: 3, 4.5, 6, 7.4, 9. At these values 
of the pH, qpip is roughly -1.5, -2.7, -3.2, -4.2, -5.0, 
respectively, based on acid dissociation constants from 
Ref. [IT]- However, because the ionization state of PIP2 
may be influenced by various geometric and chemical fac- 
tors [17], we do not assume that these qpjp values are 
exact for our system. 

Simulations — We retain only the competition between 
electrostatic interactions and excluded volume repulsions 
by adopting a model in which both lipids and small ions 
are represented as charged spheres (radius Ri) with an 
excluded volume interaction given by the purely repulsive 
(truncated at its minimum and shifted) Lennard-Joncs 
potential (the WCA potential [IB])- Parametrized by an 
energy scale e = kpT = 1 (our unit of energy) and length 
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Rj, this potential takes the following 



form as a function of ccnter-to-center distance r 
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VwcA,ij(nj) = 4e 
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for rij < 2 1 / 6 a i j, and V(rij) — otherwise. Note that cr^ 
is the distance at which the potential equals kpT. N — 
1600 lipid particles are confined to the z = plane, to 
mimic the effect of the hydrophobic interaction that keeps 
them at the air- water interface. We use Ri = R\, = 3 A 
for the lipids and Ri = Rci = 2 A for the small cations 
that are allowed to explore the entire simulation box. In 
a study of the dependence of the clustering on cation size, 
we vary it in the range 0.5 A < Rqi < 2.5 A. The box is 
periodic in x- and y-directions (size L x = 320AxL y = 
320 A) and has hard walls at z — and z = L z = 200 A. 



The typical distance between lipids in the monolayer at 
z = is therefore 8 A. 

In addition, the charged spheres interact via the 
Coulomb interaction, Vc,ij — qiqjlB/fij, where we mea- 
sure charges q in units of the proton charge. In room 
temperature water the Bjerrum length l& « 7 A. 

We run molecular dynamics simulations using 
LAMMPS Q2], with a Nose-Hoover thermostat [20] and 
PPPM for the long range Coulomb interactions [2"T] . 

The strong Coulomb attraction between the anionic 
lipids and the small cations allows them to bind at a 
distance of rougly er^-. The essence of ion-mediated at- 
tractions is that these bonds are strong and long-lived 
enough so that one or two counterions can draw together 
two lipids and be bound to both simultaneously [2*2*] . 
Due to its coarse-grained nature, our model underesti- 
mates the binding energy of such bonds in two ways. 
Firstly, in real PIP2 the negative charges are localized 
mainly in phosphate groups that lie close to the surface 
of the molecule, so that the distance between the phos- 
phate groups and the cations is much closer than the 
typical 5 A allowed by our "spherical" lipids. To deter- 
mine the effect of the inter-charge distance on the bind- 
ing energy, we performed two test calculations with two 
"lipids" (<7pip = —4) and four divalent ions. We found 
that the binding energy is 2.2 times larger for a struc- 
tured lipid with two charges (q — —2) about 4 A apart, 
than for the spherical lipids in our model. Secondly, Ca 2+ 
can be expected to lose some of its hydration shell when 
it binds to a phosphate group, so that the effect of water 
as a dielectric medium is partially eliminated. This fur- 
ther increases the net binding energies by an unknown 
amount. As an estimate, we assume the binding energies 
produced by our model interaction are a factor of 3 lower 
than those in the real system and compensate for this by 
using a dielectric constant lower by a factor of 3; i.e., we 
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FIG. 2: Phase diagram (charge vs. PlP2-fraction) of a 
charged-neutral mixed monolayer (a) in experiments, where 
the lipid charge is set by the pH, and (b) in simulations (diva- 
lent ion radius Rci = 2 A). The closed discs in the shaded co- 
existence region indicate where clustering was observed. Open 
circles mark mixed samples, and grey discs are too close to 
the boundary to determine their behavior with certainty. 



FIG. 3: Dependence of the minimum magnitude of lipid 
charge needed for cluster formation on the divalent ion radius. 
Colors and symbols are the same as in Fig. [2] (a) Experimen- 
tally (<^>pip = 0.25) cluster formation occurs more easily with 
the smaller of two tested ions (Ca 2+ ). (b) The simulation 
(</>pip = 0.05) shows the same trend with ion radius. 



use a Bjerrum length of 21 A. 

Results — At high PiP2-charge, for example at pH 7.4 
in the experiment, where qpip ss —4.2, or at qpjp = —4 
in the simulation, cluster formation is readily observed 
(Fig. [l}. The experimental image (Fig. [T^,) shows the 
larger-scale picture of bright fluorescent spots marking 
the regions of large PlP2-concentration, while the simu- 
lation image (Fig. [T]b) shows still growing clusters at a 
length scale that is 1000 times smaller, after 3.5 ps of sim- 
ulation time. As expected, the positions of the condensed 
calcium ions (red discs in Fig. [iJd) clearly indicate their 
role in binding the charged lipids (green discs) together. 

The morphology observed in the early stages of coars- 
ening in the simulations illustrates some particular fea- 
tures of ion-mediated attractions that set them apart 
from simple attractive potentials. As can be seen 
in Fig. [TJd, the PiP2-rich clusters are often irregularly 
shaped, and even string-like. This occurs because the 
attraction is the net result of strong attractions (PIP2- 
Ca 2+ ) and strong repulsions (PIP2-PIP2 and Ca 2+ -Ca 2+ ) 
that can both be several tens of kp,T. Hence, a rearrange- 
ment of the lipids in a cluster typically involves energy 
barriers that are much higher than the net attraction en- 
ergies, so that evolution towards more compact shapes 
is severely hindered kinetically. In the earliest stages of 
coarsening, most domains are string-like, because for very 
small clusters such linear arrangements have the lowest 
Coulomb energy. As the domains grow, compact shapes 
become energetically favorable but are difficult to reach 
for two reasons. First, once there is a string-like cluster, 
the electric field in its neighborhood is focused towards 
the end of the string (see Fig. [IJ:), which makes it more 
likely for the next lipid to bind at the end, thus extend- 
ing the string. Second, the energy barrier for the string 
to fold onto itself is quite high. Therefore, string-like do- 
mains can persist even in the later stages of coarsening, 
as seen for example in the bottom left of Fig. [T|3. Such 
irregular domains have been seen experimentally [BJ |S] . 



To determine the conditions under which cluster for- 
mation occurs, a grid of parameter values (0pip, gpip) 
was explored. Experimentally, we identify which images 
have more bright spots than the baseline level (which we 
found to be at most two spots per frame). The result- 
ing phase diagram is shown in Fig. [2^,, where the region 
of cluster formation is shaded. In the simulations, we 
follow the coarsening dynamics by keeping track of the 
static structure factor of the charged lipids, 
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where N is the number of PiP 2 -particles. As a function 
of k = |k|, a maximum in this function at k = fc pea k indi- 
cates that the PiP2-postions are developing structure at a 
length scale 27r/fc pea k- For the more pronounced cases of 
cluster formation (deep in the phase-separated regime), 
we followed this peak as a function of time and verified 
that it scales with time as fc pca k ~ i -1 / 3 , consistent with 
the general theory of coarsening of a binary fluid mix- 
ture |23j . Thus, even though the counterion-mediatcd 
origin of phase separation yields irregularly shaped clus- 
ters instead of circular ones, this does not seem to af- 
fect the kinetics of coarsening. In the phase diagram in 
Fig. [2)3, all parameter values (</>prp, Qpip) for which an ap- 
preciable peak appears that approaches fc pca k = in S(k) 
for long times were marked as cluster-forming (within the 
coexistence region). Both in the experiment and simula- 
tion, we found that divalent cations cause phase separa- 
tion provided the lipid charge is high enough (pH 4.5 or 
higher in experiment, qpip < —2 in simulation). Mono- 
valent cations were never seen to induce clusters. 

Larger divalent ions than Ca 2+ should mediate weaker 
attractions, because larger binding distances imply lower 
Coulomb energies. This effect should manifest itself in a 
higher charge on the PIP2 needed to obtain cluster for- 
mation with larger ions. We verified this in experiments 
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using Mg 2+ , which is known to have a larger hydrated ra- 
dius than Ca 2+ , even though precise values appear to be 
lacking [25]. We tentatively indicate this size difference 
in Fig. [3^,. In agreement with this observation and the 
theoretical expectation, the ability of divalent cations to 
drive cluster formation in our simulations also decreases 
with increasing ion size (Fig. [3]b) . A subtle effect that 
could influence the effective size of the ions is the level of 
dehydration that occurs: the more of the hydration shell 
is removed, the closer and hence stronger the binding. 
We have therefore refrained from attempting to assign 
numerical values to the radii of the cations in Fig. [3^,. 

Discussion — The phase diagram of our model com- 
pares surprisingly well with our experiments without 
any parameter optimization. The only free parameter 
is the dielectric correction factor we apply to compen- 
sate for the underestimation of the mediated binding en- 
ergy, which we set to 1/3 from the start, based on cal- 
culations of the electrostatic interactions involved. The 
introduction of this correction factor allows a major sim- 
plification, of replacing lipids by spheres of the appropri- 
ate charge. This enabled us to explore a large param- 
eter space with modest computational resources. Obvi- 
ously, this non-systematic coarse-graining approach does 
not guarantee quantitative precision, but we note that 
we obtain quantitative agreement with experimental val- 
ues and trends of the surface pressure in the lipid mono- 
layer before and after adding Ca 2+ . In both simulation 
and experiment at all parameter values studied, surface 
pressures were of order 10 to 30mN/m and dropped by 
several mN/m upon addition of Ca 2+ (data not shown). 

It should be noted that, while hydrogen bonds between 
the PiP2-molecules exist and may play a role when the 
charges are small [T2], our work strongly suggests that 
they do not play a dominant role in multivalent ion- 
induced clustering — if they did, having a higher PIP2- 
charge would make it harder to form clusters, rather than 
easier, as we report in Fig. [2j 

Since the interactions in our model have been stripped 
down to the bare minimum of electrostatics and steric 
repulsion, the only attractive interaction in the simu- 
lations is the Coulomb attraction between the anionic 
lipids and the divalent cations. Therefore, the observed 
phase separation must be due to cation-mediated attrac- 
tions. The fact that these attractions are the result of 
near-cancellation of even larger attractions and repul- 
sions leads to several special features. We have already 
mentioned the high energy barriers for lipid rearrange- 
ments that lead to long lived irregular or string-like do- 
main shapes in our simulations. Such domains are ob- 
served experimentally [6] at pH values close the physio- 
logical value (where qpip is high). It is therefore possible 
that stringlike clusters appear in biological contexts. 

Another striking aspect of ion-mediated attractions is 
the strong dependence of the effective attraction strength 
on the lipid charge. To illustrate this point, we numeri- 



cally calculated the binding energy per lipid for a cluster 
of 30 lipids [26], and found that it increases from ik^T 
at gpip = —2 to 6/cbT at gpip = —3 and 11 fee? 1 at 
<Zpip = — 5. At grpip = —2 phase separation first appears, 
and at gpip = —3 it is already strong enough to lead 
to nearly complete phase separation (see Fig. [2]). This 
has the potentially important biological consequence that 
clusters can form even at the extremely low concentra- 
tions of pip 2 (~ 1%) found in the cell membrane. 
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